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We present recent D0 results based on approximately 1 fb — 1 of pp collisions at \/s = 1.96 TeV recorded 
at the Fermilab Tevatron. Preliminary results on a search for the flavor changing neutral current process 
D + —* n + fi + ii~ , a measurement of the CP violation parameter in B mixing, eg, and a two sided limit on the 
B s oscillation frequency Am s are presented. The limits on eg and S(D+ — > 7r+/^+/i~) are the world's best 
limits. The two sided bound on Am s is the first direct indication by a single experiment that Ara s is bounded 
from above. 
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1. Introduction 

A major goal of the heavy flavor physics community 
is to find indirect evidence of new physics by study- 
ing neutral meson mixing and flavor changing neutral 
current (FCNC) decays. In the Standard Model (SM) 
these processes are mediated by higher order weak 
transitions and can receive sizeable corrections from 
new particles as predicted by several extensions to the 
SM. These new physics corrections are manifested in 
either an augmentation of the decay or mixing rate, or 
possibly by introducing new complex operators lead- 
ing to unexpected sources of CP violation. 

The D0 experiment operating at the Fermilab 
Tevatron has a rich heavy flavor program built on 
the excellent performance of the D0 muon system. 
The large single and dimuon data samples are ideal 
for studying heavy flavor mixing and radiative de- 
cay. In this presentation, we focus on three re- 
sults that highlight the strengths of the program: a 
search for the flavor changing neutral current transi- 
tion c — > ufj, + /i~ []], a search for CP violation in B 
mixing [2j , and our progress in measuring the B s mix- 
ing frequency ,3'j. These results are derived from a 
data sample of pp collisions at y/s = 1.96 TeV corre- 
sponding to a luminosity of approximately 1 fb _1 . 



2. D0 experiment 

The D0 detector is a general purpose spectrome- 
ter and calorimeter Q. Charged particles are recon- 
structed using a silicon vertex tracker and a scintil- 
lating fiber tracker located inside a superconducting 
solenoid coil that provides a 2 T magnetic field. Pho- 
tons and electrons are reconstructed using the inner 
region of a liquid argon calorimeter optimized for elec- 
tromagnetic shower detection. Jet reconstruction and 
electron identification are further augmented with the 
outer region of the calorimeter optimized for hadronic 
shower detection. Muons are reconstructed using a 
spectrometer consisting of magnetized iron toroids 
and three super-layers of proportional tubes and plas- 
tic trigger scintillators located outside the calorimeter. 



The D0 trigger is based on a three tier system. The 
level 1 and 2 dimuon triggers rely on energy deposited 
in the muon spectrometer and fast reconstruction of 
muon tracks. For single muon triggers the muon is 
required to be associated with a track in the central 
detector with a transverse momentum of at least 3 
GeV/c. The level 3 trigger performs fast reconstruc- 
tion of the entire event allowing for further muon iden- 
tification algorithms, improved matching of muon can- 
didates to tracks reconstructed in the central tracking- 
system, and requirements on the z position of the pri- 
mary vertex. At higher luminosities, we increase the 
penetration requirements on the dimuon triggers, the 
momentum thresholds are usually raised to as high as 
5 GeV/c on the single muon triggers, and invariant 
mass or displaced particle filters are added at level 3 
to enhance the heavy flavor content of the muon sam- 
ples. 



3. Flavor Changing Neutral Current 
Charm Decays 

The excellent agreement between observed FCNC 
processes involving down- type charge —1/3 quarks 
such as b — > 57, b — > sl + l~, and K — > ~kvv with 
SM expectations have been used to set strict limits on 
new phenomena |MSI3- However, there are several 
scenarios of new phenomena such as SUSY R parity 
violation in a single coupling scheme Q , or little Higgs 
models with a new up-like vector quark p| where de- 
viations from the SM would only be seen in the up 
sector. Scenarios of this nature motivate the study of 
FCNC charm meson decays. 

Due to GIM suppression, the SM rates for FCNC 
charm decays vanish in the limit of SU(3) symmetry. 
The inclusive rates of decays such as D + — ► n + l + l~ 
are therefore expected to be dominated by long dis- 
tance contributions where the dilepton system is pro- 
duced via intermediate strong resonances such as the 
4> or the u> 0, El- Of particular interest is the de- 
cay — > 7r + /i + /i~ that can only proceed via the 
long distance interaction where the dimuon system 
is produced by an intermediate (f> meson. Since this 
has no short distance contributions, the rate is given 
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simply by the product of the two branching fractions 
of .D+ — > </>7r + and — > Just as the de- 

cay £? + — > J/ij}K + — > K + l + l~ played a crucial role 
in benchmarking the studies of b — > sl + l~ transi- 
tions, the observation of the Cabbibo favored decay 
Df — > </)7T + — > is an essential first step in 

the study of c — > transitions. Once this mode 

has been observed, the search for continuum produc- 
tion of c — * ul + l~ can be performed by looking for an 
excess of D + candidates where the dilepton mass is 
inconsistent with resonance production. 

D meson candidates are formed by combining two 
well reconstructed muons with a track in the same 
jet as the dimuon system and requiring the three par- 
ticle system have an invariant mass consistent with 
a D meson. If multiple candidates exist for a sin- 
gle event, the best candidate is selected based on the 
quality of the three track vertex, the transverse mo- 
mentum of the pion, and the radial distance between 
the pion and the dimuon system. Backgrounds from 
light quarks or random combinations are reduced fur- 
ther with variables that take advantage of the long D 
meson lifetime and distinct decay topology. 
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Figure 1: The m(ir + fi + fi~) mass spectrum for the loose 
cut sample in the 0.96 < m(/i + (i~) < 1.06 GeV/c 2 bin. 
The results of binned likelihood fits to the distributions in- 
cluding contributions for , D + , and combinatoric back- 
ground are overlaid on the histogram. 



The hrst goal is to observe the resonance production 
of D s — > (fm — > Trfi + fi~ . The D meson candidates for 
events where the dimuon system has an invariant mass 
consistent with a 4> meson are shown in Fig ^ for loose 
selection requirements and in Fig [21 for requirements 
optimized for resonance production and background 
suppression. In both plots a D s signal is clearly ob- 
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Figure 2: The m(7T n ) mass spectrum for the opti- 
mized Dt sample in the 0.96 < m(fi + fi~) < 1.06 GeV/c 2 
bin. The results of binned likelihood fits to the distribu- 
tions including contributions for Df , D + , and combina- 
toric background are overlaid on the histogram. 



served and in the optimized selection the D + signal is 
also seen with a significance above background of 3.1 
sigma. From this we extract the branching fraction 



>+/0 = (1.75±0.7±0.5) x 10" 



" 6 given by the product of the D + 
ji + \x~~ branching fractions. 



which can be compared to the recent CLEO-c mea- 
surement of (2.7+i'g±0.2) x 10~ 6 or the expected value 

of 1.77x10- 
and the 

Having clearly established the resonance produc- 
tion, we turn to a search for non resonant D produc- 
tion. The D meson candidates for events passing tight 
requirements optimized for the non resonant analysis 
are shown in Fig The distribution on the left are 
again events where the dimuon system has an invari- 
ant mass consistent with a <f> meson. We expect to still 
see some excess from the D s and D + — > 4>ir + signals 
and indeed based on the yields from the above sample 
we predict 5.3 ± 1.5 events (signal plus background) 
in the D + signal region and 4.6 ± 0.7 events in the 
D s signal region which compare well to the observed 
6 events in the D + window and 3 events in the D s 
window. This indicates that even after the tighter se- 
lection we still have sensitivity at the level of the D + 
branching fraction measured above. 

The distribution in Fig 0] shows the D candidates 
where the dimuon system has an invariant mass out- 
side the 4> window. Here we see 17 events in the D + 
window that are consistent with the expected back- 
ground of 20.9 ± 3.4 events. From this, we set the 
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Figure 3: The m(n (i fiT) mass spectrum for the op- 
timized D + sample in the 0.96 < m(jj, + fi~) < 1.06 
GeV/c 2 bin. The line is a fit to the data in the side- 
band regions 1.4 < m(7r + (i + fj,~) < 1.7 GeV/c 2 and 
2.1 < m^W") < 2.4 GeV/c 2 . The £>+ and D+ signal 
yields are fixed based on the yields in the optimized Df 
sample and the relative efficiencies between the Df and 
D + cuts. 



90% confidence level upper limit 

B(D+ -» 7T + /z + M~) < 4.7 x 10~ 6 . 

This is the most stringent limit to date in the dimuon 
modes. Although this is approximately 500 times 
above the SM expected signal it is already below the 
allowed parameter space of SUSY R parity violating 
couplings. 



Figure 4: Final search for D" 1 
uum region. 



7r + /i + /i in the contin- 



At D0 , we can reach the desired accuracy by tak- 
ing advantage of the two independent magnetic spec- 
trometers: the inner solenoid and the outer muon sys- 
tem toroid. The field polarities are reversed frequently 
leading to roughly equal data sets with the four pos- 
sible combinations of field configurations. Comparing 
yields in these data sets allows us to make precise de- 
terminations of all detector induced asymmetries. 

The raw event yields for different combined toroid 
solenoid polarities are listed in Table U for our 1 fb _1 
data set. The like-sign dimuon asymmetry after re- 
moving detector effects is 



.4 



-0.0013 ±0.0012 ±0.0008. 



4. CP Violation in B Mixing 

The CP eigenstates of neutral B mesons are a linear 
combination of the flavor eigenstates 

B CP =pB Q + qB°. 

A difference in the coefficients q ^ p would lead to 
CP violation B mixing. The signature of B mixing is 
like-sign dilepton events indicating that the initially 
produced BB system decayed as either a BB or a 
BB combination. The experimental signature for CP 
violation in B mixing would be an asymmetry between 
the number of l + l + and l~l~ final states. 

Since the current experimental limits are slightly 
below the percent level, our goal is to measure the 
charge asymmetry between like-sign dimuons at the 
tenth percent level precision. 



To extract the underlying CP violating semileptonic 
asymmetry A$l from this number requires precise 
knowledge of the composition of our dimuon sam- 
ple. Besides our signal of direct semileptonic decay of 
both b hadrons in the event, an important component 
of the sample is the coincidence of a primary muon 
from semileptonic b decay and a secondary muon from 
semileptonic charm decay. Topologies like these that 
produce like-sign dimuons but are not associated with 
B mixing dilute the asymmetry. 



Table 1 Raw Results. 



Toroid*Solenoid Polarity 


-1 


+1 


N++ 


177,950 


156,183 


N-- 


176,939 


156,148 


N+ - 


1,175,547 


1,029,604 
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The most problematic component are asymmetries 
induced by hadronic interactions of kaons in the de- 
tector. Reactions like K~ + N — > Air have no K + 
analog. As the kaons travel through the detector the 
fraction of K + increases and leaving a larger number 
of n~ from K + decay in flight than /i + from K~ decay 
in flight. We measure this asymmetry directly from 
reconstructed B — ► D* [iv events in the data sample 
and subtract it from the above asymmetry. This is 
by far the largest source of systematic error in the 
analysis. 

A sensitive cross check of the sample composition is 
the extraction of the integrated B mixing probability 
X that can be extracted from the ratio of like-sign 
dimuons to opposite-sign dimuons. In our data set we 
find 

X = 0.136 ±0.001 ± 0.024. 

This can be compared to the Particle Data Group 
(PDG) HH value of x(PDG) = 0.1281 ±0.0076. With 
our knowledge of the sample composition, we can 
then extract from the raw asymmetry the semileptonic 
asymmetry 

A SL = -0.0044 ± 0.0040 ± 0.0028. 

From which we extract the CP violation in B mixing 
parameter cb 

n ^ B \~ = -0.0011 ± 0.0010 ± 0.0007. 

i + ks| 2 

This measurement gives the best sensitivity to CP vi- 
olation in B mixing to date and can be compared to 
the world average value of all previous measurements 
(0.0005 ±0.0031) [13. 



5. B s Mixing 

Measurements of flavor mixing in neutral meson 
systems have historically led to profound insights into 
physics at energy scales inaccessible to particle accel- 
erators. The discovery of mixing and CP violation 
in the neutral kaon system led Kobayashi and 
Maskawa to predict the existence of a third flavor gen- 
eration 0| • The discovery of mixing in the Bd system 
gave the first indication that the top quark is much 
heavier than the W boson H3- Finally, mixing in- 
duced CP violation measurements in the Bd system 
have precisely determined the phase of the CP vio- 
lating top quark coupling Arg(V t( z) ~ 24° ^5|- This 
motivates the study of mixing in the B s system. A 
B s mixing measurement allows a precise determina- 
tion of the magnitude of the CP violating top cou- 
pling through the ratio |Vt<j/Vf S | an d may lead to the 
discovery of new physics in b — ► s transitions |16| . 



Measuring the B s oscillation frequency requires 
knowledge of the B s flavor at both production and de- 
cay. Signatures of mixing are candidates where the fla- 
vor at decay is different than the flavor at production. 
The mixing frequency can be determined by looking 
for an oscillation in the fraction of mixed events as a 
function of the B s decay length. 

The B s flavor at decay is determined by recon- 
structing a flavor-specific final state. At D0 , we use 
the semileptonic decay B® — > D~fi + v, D~ — > 4>ir~ , 
<f> — > K + K~ and charge conjugate states. The <jm in- 
variant mass distribution is shown in Fig [S] In our 
1 fb _1 sample, we reconstruct approximately 26.7 
thousand D s /i candidates predominantly from B s de- 
cay and 7.4 thousand D + fi candidates predominantly 
from B° decay. The B s flavor at production is deter- 
mined by examining inclusive properties of the sys- 
tem recoiling against the B s candidate. In this way, 
we determine the flavor of the other B hadron in the 
event that was pair produced with the B s candidate. 
We use a likelihood based tagging algorithm based on 
momentum weighted sum of the charge of particles 
associated with either a lepton, a displaced vertex, 
or the entire recoil system. Fig [SJd shows the event 
sample when the flavor tag is required. The tagging 
purity of the combined likelihood is determined event 
by event and is calibrated by studying Bd mixing in 
our B® — > D*~ sample. The figure of merit for 
tagging is the effective efficiency, eD 2 , that takes into 
account both the efficiency of the tagging algorithm 
and the weight that each event receives in the final 
analysis due to the purity of the tag. In our D* jiv 
sample, we measure an effective efficiency of 

eD 2 = 2.48 ±0.21±°;°s. 

The decay length of the B s meson is taken as the 
transverse distance from the primary vertex recon- 
structed using tracks originating in the beam spot, 
and the secondary vertex reconstructed from the D s 
and fi particles. To boost the decay length to the 
proper frame, we use MC to simulate the missing en- 
ergy from the neutrino or other pions or kaons from 
B s — ► D S X decays. 

We require a very accurate determination of our 
decay length resolution so that we can distinguish a 
lack of a signal in our data from a lack of the ability to 
resolve a frequency beyond our resolution. We test our 
resolution model using prompt J/ip candidates from 
pp — » J/ipX interactions. Here the J/ij) decays at 
the primary vertex so that our measurement of the 
distance from the primary vertex to the J / ijj vertex is 
a direct measurement of our resolution. 

Fig shows the value of —AlogC as a function of 
Ato s , indicating a favored value of 19ps , while vari- 
ations of -logC from the minimum indicates an oscilla- 
tion frequency of 17 < Am s < 21ps~ 1 at the 90% C.L. 
The uncertainties are approximately Gaussian inside 
the interval. Using 1000 parameterized Monte Carlo 
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samples, it was determined that for a true value of 
Am s = 19ps _1 , the probability was 15% for measur- 
ing a value in the range 17 < Am s < 21ps _1 with a 
—AlogC lower by at least 1.9 than the corresponding 
value at Am s = 25ps~ 1 . 

The amplitude method 33 was a l so used and the 
results are shown in Figd The unbinned likelihood fit 
is repeated including the oscillation amplitude A as a 
function of the input value of Am s . At Am s = 19ps _1 
the measured data point deviates from the hypothe- 
sis A — (A — 1) by 2.5 (1.6) standard deviations, 
corresponding to a two-sided C.L. of 1% (10%), and 
is in agreement with the likelihood results. This is 
the first report of a direct two-sided bound on the B° s 
oscillation frequency. 
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Figure 7: B® oscillation amplitude as a function of oscil- 
lation frequency, Am s . 




1.8 1.9 2.0 1.8 1.9 2.0 



Figure 5: (K + K~ ')ir invariant mass distribution for (a) 
the untagged sample, and (b) for candidates that have 
been flavor-tagged. The left and right peaks correspond 
to fi + D~ and u + D~ candidates, respectively. The curve 
is a result of fitting a signal plus background model to 
the data. For fitting the mass spectra, a single Gaussian 
function was used to describe the D~ — > 4>tv~ decays and 
a double Gaussian was used for the — > <^7r~~ decays. 
The background was modeled by an exponential function. 



6. Conclusions 

D0 preliminary results based on approximately 1 
fb -1 of pp collisions at yfs — 1.96 TeV recorded at the 
Fermilab Tevatron have been presented. We set the 
most stringent limit to date in a search for the flavor 
changing neutral current in the dimuon modes. At 
the 90% confidence level upper limit we find B(D + — * 
tt+^ + aO < 4.7 x 10~ 6 . We have extracted the CP 
violation in B mixing parameter es 

- K ^ B \ = -0.0011 ± 0.0010 ± 0.0007, 

which is the best sensitivity to CP violation in B mix- 
ing to date. We have presented the first report of a di- 
rect two-sided bound on the B° s oscillation frequency, 
17 < Am s < 21PS" 1 at the 90% C.L. 
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Figure 6: Value of —AlogC as a function of Am,. Star 
symbols do not include systematic uncertainties, and the 
shaded band represents the envelope of all logC scan curves 
due to different systematic uncertainties. 
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